Growth Models and Models of Turbulence : A Stochastic Quantization Perspective 
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We consider a class of growth models and models of turbulence based on the randomly stirred 
fluid. The similarity between the predictions of these models, noted a decade earlier, is understood 
on the basis of a stochastic quantization scheme. 
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Over the last decade, it has been noticed that two dif- 
ferent problems in non equilibrium statistical mechanics 
seem to have the same characteristics. One is the growth 
of interfaces under random deposition - the sort of growth 
that characterises molecular beam epitaxy and the other 
is the behaviour of velocity in a randomly stirred fluid, 
which arguably is a decent model of turbulence. The 
surface deposition model has its origin in the works of 
Wolf and VillainQ and Das Sarma and TamboreneaQ- 
These were discrete models where particles were dropped 
at random on an underlying lattice. Having reached the 
surface, the particle diffuses along the surface until it 
finds an adequate number of nearest neighbours and then 
settles down. The surface is characterised by the corre- 
lation of height fluctuations. The equal time correlation 
of the height variable h{r, t) is expressed as 



< [A/i(r, >=< [h{f+ ro, t) - /i(ro, t)f >(x r^" 



(1) 



as t becomes longer than L^, where L is the system 
size. In the above, a is called the roughness exponent 
and z the dynamic exponent. The surface is rough if a is 
positive. The continuum version of this model was intro- 
duced by Lai and Das SarmaQ- The similarity of this 
growth model with models of fluid turbulence was first 
noticed by KrugQ who studied the higher order correla- 
tion functions and found 



< [Ah{f)YP >cx r 



(2) 



where ap is not a linear function of p. This is a man- 
ifestation of intcrmittency - the single most important 
feature of fluid turbulencej^ IE 0- It says that the 
distribution function for | Ah{f,t) \ has a tail which 
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is more long lived than that of a Gaussian. The rele- 
vant quantity in fluid turbulence is the velocity increment 
Ava = Va{r + ro, t) — Vaifo^ t), where we are focussing on 
the a*'' component. It is well known that Ava exhibits 
multifractal characteristics or intermittency and the mul- 
tifractal exponents have been measured quite carefully 
The multifractal nature of the velocity field is expressed 
by an equation very similar to Eq.Q, 



< 



(3) 



where Cp is not a linear function of p. From some con- 
siderations of singular structures, She and Leveaue[lo| 
arrived at the formula (p = ^ + 2[1 - (1)^^/^], which 
gives a very reasonable account of the experimentally de- 
termined multifractal indices. It is the similarity between 
the growth model and the turbulence characteristics that 
is interesting. 

The randomly stirred Navier Stokes equation^J 
113 . ITU in momentum space reads 



p 



-p) =ikaP-iyk^Va + fa (4) 



with the incompressibility condition appearing as 
0. The random force fa has the equal time 
correlation 



ka'^a {k) 



< fa{k^)fli{k2) >= 2 



Do 



5ali5[ki + k2 



(5) 



where D is the dimensionality of space and j/ is a pa- 
rameter supposed to give the fully developed turbulence 
for y = 4. In this model, the pressure term does not 
qualitatively alter the physics of turbulence and hence it 
is often useful to study the pressure free model. This was 
extensively done by Checklov and Yakhot|l6j and Hayot 
and ,TavaDrakash|17j in D = 1. If we write Va — dah, 
then 



h{k) = -vk^h{k) - l/2^p- (fc- 



-p)h{p)h{k^p)+9{k) 
(6) 
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where 



< 9c{k)g/3{k') >-- 



2Do 



2D, 



kD-2+ySik + k') = ^6{k + k') (7) 

This is the Medina, Hwa and Kardar model[l3 in 
D-dimensions. Turning to the growth model, the Un- 
ear equation for surface diffusion hmited growth is the 
Mulhns - Sereska modellS*] given by 



dhjk) 



(8) 



where, < rj{k)rj{k') >= 2DQ6{k + k') and its generali- 
sation to include nonlinear effects is the Lai-Das Sarma 
niodel0| defined as 



dh{k) 
dt 



-yk%{k)- 



k'^^p-{k-p}h{p)h{k-p)+r](k) 



(9) 

The various properties of the model have been very 
well studiedjM M IH El Hi • 

Our focus is on the similarity between the growth 
model and the model of turbulence. In some sence, these 
two widely different models (one with coloured noise, the 
other with white noise ) have to be related. We introduce 
a technique of handling non equilibrium problems that is 
based on stochastic quantization|23 and show that the 
two models can be made to look quite similar. We con- 
sider a scalar field 0(fc) in momentum space satisfying 
the equation of motion 



m^-mm-M{4>)+9{k) 



(10) 



M{(j)) is a non linear term in (/) and g{k) is a noise 
which can be coloured in general and we will take it to 
be of the form of Eq.©. The probability distribution 
corresponding to the noise is given by 

f d^k dw k^P , , , N , ^ 

Piqjoaexp— / - — -=r q{k,w)q{—k,—w) (11) 

^^'^ ^ J {2t:)° 27r AD^^ ' , I \ ) 

In momentum- frequency space, Ea. (|10|) reads 

[-IW + i(/c)]0(fc, w) + Mfc,^ (0) ^g{k,w) (12) 

The probability distribution written in terms of (/)(fc, w) 
instead of g{k, w) is 



P oc exp{ 



d^k dw 



ADo J (27r)^ 2tt 
{[-iw + L{k)]^{k, w) + Mk,M)} 
{[iw + L{-k)]cj){-k, -w) + M_k,-M)}} 



V[4>]e 



-J— f 



'-S(k,w) 



(13) 



At this point of development, the usual practice is to 
introduce a response field (f>, work out the response func- 
tion 'AS < 4>(p > and the correlation function as < >. 
There is no fiuctuation dissipation theorem to relate the 
two and hence two independent considerations are neces- 
sary. We now exploit the stochastic quantization scheme 
of Parisi and Wu to introduce a fictitious time V and 
consider all variables to be functions of r in addition to 
k and w. A Langevin equation in V space as 



d4>{k, w, t) 



5S 



'k, — u), t) 



+ n{k,w,T) (14) 



with < nn >= 2S{k + k')S{w + w')S{t - r'). 

This ensures that as r ^ cxd, the distribution function 
will be given by S{k,w) of Ea. (|13|l and in the r-space 
ensures a fluctuation dissipation theorem. From Ea. ((T^ . 
we find the form of Langevin equation to be 



d4>{k, w, ' 



—p.—w' (</>) 



d^p dw' 



p'^ii-i', 

+ {iw' + L{-p))(j){-p,w')Mj;^^,]] 

5 , f d^p dw' 



(27r)^ 2tt 



■p"'M^^^,{cp)M_^,. 



5(t)'J (27r)^ 27r 

'((/))] -|-n(fc,w,T) 
(15) 



The correlation functions calculated from the above 
Langevin equation lead to the correlation functions of 
the original model as r ^ oo. For proving scaling laws 
and noting equivalences, it suffices to work at arbitrary 
T. It is obvious from Ea. ljlSI) that in the absence of the 
nonlinear terms (the terms involving M {(/))), the Greens 
function G*^°^ is given by 



[G^"y = -inr + k^''- 



2Dn 



(16) 



where fir is the frequency corresponding to the fictitious 
time r. As is usual, the effect of the nonlinear terms, 
leads to the Dysons' equation 



G-i = [GW]-i+^(fc,«;,n.) 



(17) 



The correlation function is given by the ffuctuation dis- 



sipation theorem as C 



-IraG. Let us start with 



Eq.@ which relates to ffuid turbulence. The linear part 
of the corresponding Eq. Hlt)|) gives 

V.2p 

[G(°)]-i = -zl7. + |-(^2+z.2fc4) (18) 



The 



Qo limit of the equal time correlation func- 



tion is 2Do/k'^''{w'^ + f^k'^), which leads, in D = 1, to 
a = {1 + 2p)/2. The dynamic exponent is clearly z = 2. 
If we now turn to the growth model of Eq.® and con- 
sider the linear part of the relevant form of Ea. (|15|l . then 
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[-(^(o)]-i ^ _|_ (-^2 _|_ i,2i^8y2D^ and the correspond- 

ing a = 3/2 in Z? = 1. The dynamic exponent z is 4. We 
note that although the dynamic exponents never match, 
the two roughness exponents are equal for p = I. This is 
what is significant. 

We now turn to the nonlinear terms and treat them to 
one loop order. For Eq.©, Mk,wi(l>) = (1/2) EpP' " 
p)(bip)(j){k-p), while for Eq.® Mk,M) = {1/2)P ^-p- 

{k — p)(j){p)(j){k — p). The nonlinear term which involves 
two Af's in Ea. (fTK|l gives a one loop correction which 
is independent of external momenta and frequency and 
hence is not relevent at this order. It is the term involv- 
ing one M which is important and for Eq.®, this has 
the structure k'^P {—iw + vk^){p- [k — ■p))(j){p)(j)(k —p). For 
Ea. (|10|l . the corresponding structure is {—iw + i'k'^)[k'^p- 
[k — p)](j){p)(j){k — p). For p — 1, the two nonlinear terms 
have very similar structure! The scaling of the correla- 
tion function determines the roughness exponent. Now 
the one loop graph in both cases are composed of two ver- 
tices, one response function and one correlation function. 
While the dynamic exponent z will differ the momentum 
count of the one loop graph for the fluid must agree with 
that for the interface growth since for p — 1, the ver- 
tex factor agree, the correlation functions tally and the 



frequency integrals of G*^"^ match. Thus at p = 1, the 
perturbation theoretic evaluation of a for for the two 
models will be equal. 

How big is a in the growth model? A one loop self 
consistent calculation yields the answer in a trivial fash- 
ion. The structure of is / d'^pdwdVtrV.V.GC . We 
recall that C is \/VlrIrnG and hence dimensionally this 
is / d^pdwVVGG. If the frequency scale is to be modi- 
fied to fc^ from fc^ with z < 4, then G scales as and 
V fc'*+^ and hence ^ k^^^~^ which has to match k"^^ . 
This yields z = {D + 8)/3. 

A ward identity shows a + z — A and thus a = 
(4 ~ D)/?,. At D = 1, a = 1 and matches the p = 1 
results of a = 1 for the fiuid model. But as is apparent 
from the work of referrence(17) , this is where the mul- 
tifractal nature sets in for the fluid because of the non- 
linear term. The identical structure of the growth model 
nonlinearity tells us that in D = 1, it too will have mul- 
tifractal behaviour. Thus, we see that the growth model 
and the turbulence model are not in the same universality 
class since the dynamic exponents are different but the 
structure of the Langevin equation in the fictitious time 
makes it clear that they will have the same roughness 
behaviour. 
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